AD-A103  417  ELECTRONICS  RESEARCH  LAB  ADELAIDE  (AUSTRALIA)  F/S  20/14 

IONOSPHERIC  TRANSEOUATORIAL  RAOIO  WAVE  OU I DANCE  AT  VHF.(U) 

JAN  SI  C  WINKLER 

UNCLASSIFIED  ERL-0181-TR  _  NL 


AD  A 1 0  3  4 1 7 


ERL0181-TR  V 


L:  1-1  ■ 


it*  ■ 1 


/ 


AR-002S41 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
ELECTRONICS  RESEARCH  LABORATORY 

DEFENCE  RESEARCH  CENTRE  SALISBURY 
SOUTH  AUSTRALIA 


) 


TECHNICAL  REP«T  .  a 

!  ERL-0181-TR 

*  l 

/  - 

IONOSPHERIC  TR ANSEQU ATORI AL  RADIO  YVAVE  GUIDANCE  AT  VHF  , 


C.  WINKLER 


THE  UNITED  STATES  NATIONAL 
TECHNICAL  INFORMATION  SERVICE 
IS  AUTHORISED  TO 
REPRODUCE  AND  SELL.  THIS  REPORT 
I _ 


DTIC 


D 


4  ** 

±  J 


Approved  for  Public  Release 

0 

)  Commonwealth  of  Australia 

JANUAMT  t881 

80  8  28  ( 

f2%~ 

"”TJ 


The  official  documents  produced  by  the  Laboratories  of  the  Defence  Research  Centre  Salisbury 
are  issued  in  one  of  five  categories:  Reports,  Technical  Reports,  Technical  Memoranda,  Manuals  and 
Specifications.  The  purpose  of  the  latter  two  categories  is  self-evident,  with  the  other  three  categories 
being  used  for  the  following  purposes: 


Reports 


documents  prepared  for  managerial  purposes. 


Technical 

Reports 


records  of  scientific  and  technical  work  of  a  permanent  value  intended  for  other 
scientists  and  technologists  working  in  the  field. 


Technical  intended  primarily  for  disseminating  information  within  the  DSTO.  They  are 

Memoranda  usually  tentative  in  nature  and  reflect  the  personal  views  of  the  author. 


Acc'ww?cht  7xrr 

Nllir  (VAil 
FAR 

Uhawnorajcstf  n- 

i  coA ’on. _  _  _ _ _ _ 


9r _ 

V-i  ^r^or/* 
j . r  r c. it y  0 »>•!♦» 
1  ’  i  iiz'.i/u r 
Dist  Kpoc ial 


1 


a 


UNCLASSIFIED 

DEPARTMENT  OF  DEFENCE 
E FENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
ELECTRONICS  RESEARCH  LABORATORY 

TECHNICAL  REPORT 
ERL-0181-TR 


AR-002-S41 


IONOSPHERIC  TRANSEQUATORIAL  RADIO  WAVE  GUIDANCE  AT  VHF 

C.  Winkler 


SUMMARY 

The  mechanism  of  reflection  of  radio  waves  at  VHF  from 
irregularities  aligned  with  the  earth's  magnetic  field  in 
the  equatorial  region  is  examined  in  detail  and  a  computer- 
based  study  of  the  resulting  phenomenon  is  presented. 


DTIC 


SELECTl 
AU3  28  1981 


D 


POSTAL  ADDRESS:  Chief  Superintendent,  Electronics  Research  Laboratory, 
Box  2151,  GPO,  Adelaide,  South  Australia,  5001. 


UNCLASSIFIED 


DOCUMENT  CONTROL  DATA  SHEET 


Security  classification  of  this  page 


DOCUMENT  NUMBERS 


AR 

Number:  AR-002-541 


UNCLASSIFIED 


SECURITY  CLASSIFICATION 


ERL-0181-TR 


Report 

Number: 


Other 

Numbers: 


3  TITLE 

IONOSPHERIC  TRANSEQUATORIAL  RADIOWAVE  GUIDANCE  AT  VHF 


a.  Complete 
Document : 

Unclassified 

b.  Title  in 

Unclassified 

Isolation: 

c.  Summary  in 
Isolation: 

Unclassified 

PERSONAL  AUTHOR(S): 


C.  Winkler 


7.1  CORPORATE  AUTHOR(S): 


Electronics  Research  Laboratory 


7.2  DOCUMENT  SERIES 
AND  NUMBER 

Electronics  Research  Laboratory 
0181-TR 


10  IMPRINT  (Publishing  organisation) 


Defence  Research  Centre  Salisbury 


RELEASE  LIMITATIONS  (of  the  document): 


Approved  for  public  release 


5  DOCUMENT  DATE: 


January  1981 


6 

6.1  TOTAL  NUMBER 

28 

OF  PAGES 

6.2  NUMBER  OF 

16 

REFERENCES: 

8  REFERENCE  NUMBERS 


a.  Task:  DST  74/023 


b.  Sponsoring 
Agency: 


9  I  COST  CODE: 


174  844 


1 1  COMPUTER  PROGRAM(S) 
(Title(s)  and  language(s)) 


12.0  OVERSEAS  NO 


Security  classification  of  this  page: 


i 


Security  classification  of  this  page: 


UNCLASSIFIED 


1 3  ANNOUNCEMENT  LIMITATIONS  (of  the  information  on  these  pages): 


14  DESCRIPTORS: 

a.  EJC  Thesaurus 
Terms 


b.  Non-Thesaurus 
Terms 


Equatorial  zone 
Very  high  frequencies 
Telecommunication 

Radio  waves  Computerized  simulation 

Ionospheric  propagation 
Mathematical  models 


Transequatorial  propagation 
Geometric  opties 


lb  LIBRARY  LOCATION  CODES  (for  libraries  listed  in  the  distribution): 


17J  SUMMARY  OR  ABSTRACT: 

(if  this  is  security  classified,  the  announcement  of  this  report  will  be  similarly  classified) 


The  mechanism  of  reflection  of  radio  waves  at  VHF  from  irregularities  aligned 
with  the  earth's  magnetic  field  in  the  equatorial  region  is  examined  in  detail 
and  a  computer-based  study  of  the  resulting  phenomenon  is  presented. 


Security  classification  of  this  page: 


UNCLASSIFIED 


ERL-0181-TR 


TABLE  OF  CONTENTS 

1.  INTRODUCTION 

2 .  THE  MODELS 

2.1  The  Ionospheric  Models 

2.2  The  Ray  Approximations 

2.2.1  Ray  concepts 

2.2.2  Ray  tube  concepts 

2.2.3  A  model  for  attenuation 

2.3  A  model  for  Received  Power 

3.  RESULTS 

4.  CONCLUSIONS 
NOTATION 
REFERENCES 

LIST  OF  APPENDICES 

I  STRUCTURE  OF  THE  EQUATORIAL  FIELD 
II  MODEL  FOR  ELECTRON  DENSITY 
III  CHANGE  IN  REFRACTIVE  INDEX  AT  AN  IRREGULARITY 
IV  ANTENNA  POLAR  DIAGRAM 

V  THE  GEOMETRY  OF  REFLECTION  FROM  AN  ELLIPTICAL  SURFACE 

TABLE  1.  COORDINATES  OF  DARWIN  AND  YAMAGAWA 

LIST  OF  FIGURES 

1.  Power  from  a  transmitter  at  A  in  one  hemisphere  is  received  at  B  on 
the  opposite  side  of  the  geomagnetic  equator 

2.  Sketch  illustrating  the  formation  of  flux  tubes  from  rays  drawn  through 
adjacent  points  on  a  wavefront  in  (a)  a  graded  index  and  (b)  a  step 
index  dielectric  waveguide 

3.  Illustration  of  the  construction  of  flux  tubes  in  (a)  free  space, 

(b)  where  guidance  is  involved 

4.  Plot  of  normalised  received  signal  strength  against  distance  from  the 

equator  for  X  =  3.0  metres  and  h  =  600  kilometres 


Page  No. 
1 
2 
2 
3 
3 

3 

4 
6 
8 
9 

12 

10 


ERL-0181-TR 


5.  Plot  of  normalised  received  signal  strength  against  distance  from  the 

equator  for  X  =  4.0  metres  and  h  ■  600  kilometres 

6.  Plot  of  normalised  received  signal  strength  against  distance  from  the 

equator  for  X  =  4.5  metres  and  h  =  600  kilometres 

7.  Plot  of  normalised  received  signal  strength  against  distance  from  the 

equator  for  X  =  5.0  metres  and  h  =  600  kilometres 

8.  Plot  of  normalised  received  signal  strength  against  distance  from  the 

equator  for  X  =  6.0  metres  and  h  =  600  kilometres 

9.  Plot  of  normalised  received  signal  strength  against  frequency  for  both 
theoretical  and  experimental  results  as  discussed  in  the  text 

10.  Geometry  for  fitting  an  ellipse  to  the  dipole  field  line 

11.  Plot  showing  relation  between  electron  density  and  ionospheric  height 
for  an  u  chapman  layer  using  various  values  of  scale  height  H 

12.  Geometry  for  Antenna  Polar  Diagram 

13.  Geometry  for  reflection  of  a  ray  from  an  elliptic  surface 


1 


ERL-0181-TR 


1.  INTRODUCTION 

There  has  been  a  continuing  interest  in  the  use  of  the  VHF  band  to  provide 
communication  channels  which  supplement  those  of  the  more  conventional  but 
congested  HF  band.  The  availability  of  these  supplementary  channels  for  large 
range  transequatorial  circuits  was  first  reported  by  radio  amateurs  and  studies 
of  the  phenomenon  were  largely  experimental.  The  physical  mechanisms  involved 
in  producing  these  channels  have  been  and  still  are  poorly  understood.  It  is 
the  purpose  of  this  report  to  lay  down  the  elements  of  a  simple  theory  to 
explain  the  large  signal  strengths  observed  on  these  circuits  and  to  describe 
the  results  of  a  computer  model  used  to  calculate  them.. 

A  qualitative  review  of  a  variety  of  proposed  mechanisms  involved  in  producing 
these  long  distance  circuits  appeared  in  1969(ref.l).  This  review  also  contains 
an  excellent  historical  background  and  we  will  not  attempt  to  repeat  that  here. 
There  were  essentially  four  possible  ways  in  which  these  circuits  could  be  set 
up. 

(1)  Refraction  in  an  undisturbed  F  region  plasma.  The  occurrence  of  bulges 
in  the  ionization  densities  on  each  side  of  the  equator  is  assumed  to  provide 
sufficient  gradation  in  the  refractive  index  to  allow  this  to  take  place. 

(2)  Guidance  by  steep  field  aligned  gradients  in  electron  density.  The  ray 
path  is  assumed  to  be  given  a  curvature  equal  to  the  field  curvature  so  that 
guidance  occurs  along  steeply  varying  gradients  of  ionization. 

(3)  Scatter  from  F  region  irregularities.  An  ensemble  of  scatterers  at 
ionospheric  height  cannot  return  sufficient  power  to  the  opposite  hemisphere 
unless  it  is  shaped  to  provide  a  degree  of  coherence  that  is  deemed  to  be 
unlikely. 

(4)  Double  scattering  between  irregularities.  Scattering  from  irregularities 
on  one  side  of  the  equator  to  a  similar  set  in  the  opposite  hemisphere  provides 
the  correct  geometry  but  considerable  debate  exists  over  the  magnitude  of  the 
received  power. 

None  of  these  mechanisms  was  entirely  adequate  to  explain  the  phenomena  of 
transequatorial  propagation  (TEP)  and  since  the  study  was  of  a  semi-empirical 
nature  no  firm  conclusions  could  be  reached.  An  attempt  has  since  been  carried 
out  to  determine  the  physical  extent  of  field  guided  F  region  irregularities  at 
the  equatorial  peak  by  illuminating  them  with  a  sideways  looking  VHF  radar (ref. 2) 
but  the  nature  of  the  irregularities  remained  obscure. 

More  recently  large  scale  regions  of  depleted  ionospheric  plasma,  called  bubbles, 
have  been  discovered  and  observed  from  the  ground(ref . 3)  and  from  satellite 
(ref. 4).  A  description  of  TEP  of  radio  waves  at  VHF  has  also  appeared 
attributing  it  to  guidance(ref . 5)  along  field  aligned  tubes  of  depleted 
ionization  resulting  from  the  post  sunset  creation  of  equatorial  plasma  bubbles. 

A  large  ensemble  of  these  bubbles  were  then  imagined  to  contribute  the  observed 
characteristics  but  some  unusual  and  dubious  properties  were  ascribed  to  the 
waveguide  tubes. 

In  this  report  we  consider  the  properties  of  a  curved  idealized  dielectric 
discontinuity  formed  at  the  upper  edge  of  a  rising  plasma  bubble.  The  equations 
for  reflection  and  guidance  of  electromagnetic  waves  along  the  discontinuity  are 
programmed  onto  a  high  speed  digital  computer  and  profiles  for  the  received  power 
at  points  on  the  surface  of  the  earth  in  the  receiving  hemisphere  are  derived. 

A  comparison  with  some  experimental  results  for  a  Darwin  to  Yamagawa  circuit  is 
also  included. 


2 .  THE  MODELS 


The  situation  is  illustrated  in  figure  1  which  shows  a  transmitter  at  point  A  in 
one  hemisphere  and  a  receiver  at  point  B  in  the  opposite  hemisphere.  The 
geographic  equator  is  located  at  G  and  the  geomagnetic  equator  is  at  position  M. 
The  magnetic  field  structure  is  assumed  to  be  identical  with  that  of  a  magnetic 
dipole  located  at  the  centre  of  the  earth  so  that  field  lines  are  symmetrical 
about  the  geomagnetic  equator. 

In  the  situation  we  have  here  a  dielectric  discontinuity  is  assumed  to  be  exactly 
aligned  with  the  particular  field  line  coincident  with  the  top  edge  of  a  rising 
plasma  bubble.  The  plasma  bubble  represents  a  region  of  depleted  ionization 
(refs. 3,4,6)  for  which  we  take  the  refractive  index  as  being  identically  one. 
Where  the  ionization  is  not  totally  depleted  then  the  refractive  index  is  only 
close  to  one  but  detailed  studies  of  this  marginal  effect  are  left  to  a  later 
report. 

Since  the  complete  modal  solution  for  this  problem  is  extremely  difficult,  if  not 
impossible,  then  it  is  necessary  to  resort  to  simplified  methods.  One  method 
that  we  can  use  in  this  situation  is  based  on  some  simple  concepts  from  geometric 
optics.  We  discuss  the  ray  approximations  in  some  detail  in  the  following 
sections. 

2.1  The  Ionospheric  Models 

It  is  shown  in  Appendix  I  that  the  field  lines  near  the  equator  assume  a  shape 
very  close  to  elliptical  (within  1%).  Furthermore  it  has  also  been  shown 
that  elliptical  multimode  dielectric  waveguiding  structures (refs . 7 ,8)  can  be 
well  represented  using  an  extended  form  of  geometric  optics  which  allows  for 
tunnelling  from  curved  dielectric  surfaces. 

An  appropriate  model  for  the  ionospheric  electron  density  is  the  Chapman  model 
described  in  Appendix  II.  The  change  in  electron  density  is  related  to  a 
change  in  refractive  index  as  given  in  Appendix  III.  Combing  these  two 
results  one  can  readily  see  that  the  refractive  index  discontinuity  varies 
with  height  above  the  earth's  surface.  The  discontinuity  gives  rise  to  a 

critical  angle  0^  for  rays  incident  from  the  more  dense  medium.  We  define 


0  =  arcos  —  (1) 

c  n:  J 

where  nt  is  the  refractive  index  of  the  more  dense  medium  (=1  in  this  case) 
and  n:  is  the  refractive  index  of  the  ionized  medium. 

Classical  geometric  optics  predicts  that  rays  incident  at  less  than  0^  are 

totally  internally  reflected  while  rays  incident  at  greater  than  6  are 

refracted  according  to  Snell's  and  Fresnel's  Laws.  Since  the  electron 
density  varies  with  height,  then  the  critical  angle  is  also  height  dependent 
and  varies  along  the  magnetic  flux  line  which  contains  the  discontinuity. 

However,  in  this  case  because  of  the  curvature  of  the  magnetic  field  line  we 
must  extend  the  concepts  of  geometric  optics  to  include  the  phenomenon  of 
electromagnetic  tunnelling.  We  discuss  this  in  detail  in  subsequent  sections. 


At  heights  below  100  km  the  discontinuity  is  considered  to  be  so  small  as  to 
be  negligible. 
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2.2  The  Ray  approximations 
2.2.1  Ray  concepts 

An  electromagnetic  wave  propagating  in  free  space  has  a  well  defined 
direction  of  propagation  specified  by  wave  vector  k  normal  to  the  electric 
and  magnetic  field  vectors,  E  and  H  respectively.  The  vectors  E,  H  and  k 
form  a  right-handed  triad.  In  a  uniform  homogeneous  medium,  mT’dif ficulty 
exists  in  the  specification  of  the  direction  of  the  wave  and  the  concept  of 
a  ray  is  easily  managed.  A  ray  indicates  the  direction  of  the  local  plane 
wave  vector  k,  which  is  also  the  direction  of  power  propagation.  Maxwell's 
equations  are  exactly  satisfied  and  the  direction  of  the  wave  is  simply  the 
direction  of  a  ray. 

In  an  inhomogeneous  medium,  where  the  electromagnetic  fields  are  not  in 
general  specified  by  a  single  plane  wave,  the  direction  of  a  ray  is  no 
longer  so  simply  identified  and  it  is  necessary  to  specify  it  asymptotically 
as  the  short  wavelength  (X  +  0)  approximation  for  the  local  plane  wave 
direction.  It  is  used  to  indicate  the  direction  of  k  and  the  magnitude, 
given  by  the  relation 


Ikl 


k  = 


2irn 
X  * 


(2) 


is  taken  to  the  limit  k~*  ~.  X  is  the  wavelength  of  light  in  vacuum.  The 
ray  path  is  then  always  defined  by  the  eikonal  equation. 


V  n 


(3) 


where  s  is  the  distance  measured  along  the  path  and  is  the  position 
vector  for  a  point  on  the  ray  path.  n  is  the  refractive  index  of  the 
medium  at  the  point  under  consideration  and  d$?/ds  is  interpreted  physically 
as  a  unit  vector  tangent  to  the  ray  path.  This  is  the  asymptotic  theory 
known  as  geometric  optics.  The  ray  direction  defined  in  the  short  wave¬ 
length  approximation  is  specified  everywhere  including  caustics  and  points 
of  reflection. 

Indiscriminate  use  of  geometric  optics  leads  to  strange  contradictions.  It 
predicts  the  field  intensity  at  a  focus  or  caustic  to  be  infinite  which  is 
obviously  ridiculous.  The  field  distribution  near  caustics  is  not 
properly  calculated  using  geometric  optics  and  in  order  to  allow  the 
advantages  of  the  simplicity  of  geometric  optics  to  be  exploited  everywhere 
we  aim  to  extend  geometric  optics  to  accommodate  these  effects.  A  similar 
approach  has  already  been  used  by  Keller(ref . 9)  in  the  Geometric  Theory  of 
Diffraction  which  has  proved  highly  successful. 

2.2.2  Ray  tube  concepts 

Consider  the  situation  shown  in  figure  2  where  two  rays  are  drawn  through 
adjacent  points  on  a  wavefront  within  (a)  a  graded  index  and  (b)  a  step 
index  dielectric  waveguide  forming  a  characteristic  flux  tube  as  shown  in 
each  case.  The  flux  or  ray  tube  can  be  traced  as  shown  through  an  optical 
medium  by  tracing  each  of  the  rays  individually. 
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The  region  of  depleted  ionization  we  can  regard  as  being  the  core  of  a 
dielectric  waveguide  while  the  region  above  it  is  exterior  to  the  guide  and 
has  commonly  been  referred  to  as  the  cladding.  When  the  core  is  non 
absorbing,  power  in  the  ray  tube  travels  parallel  to  the  ray  direction  and 
can  only  be  lost  to  the  cladding  at  reflections  from  the  core  cladding 
interface  or  at  caustics  where  the  asymptotic  description  of  geometric 
optics  is  inadequate. 

2.2.3  A  model  for  attenuation 

The  definition  of  a  ray  allows  us  to  associate  with  a  ray  power  travelling 
in  the  same  direction.  Strictly  speaking,  the  terminology  should  refer  to 
either  a  small  ray  tube  or  flux  tube  but  it  is  simpler,  physically 
appealing  and  unambiguous  to  merely  refer  to  it  as  the  power  of  the  ray. 

In  a  guiding  structure  the  local  plane  wave  fields  of  a  family  of  rays,  all 
with  the  same  ray  invariants  and  following  similar  ray  paths,  is  the 
geometric  optics  equivalent  of  a  mode (ref. 7).  In  order  to  develop  a  ray 
model  for  attenuation  it  is  convenient  to  introduce  the  concept  of  a  ray 
transmission  coefficient.  If  the  refractive  index  of  the  core  is  nj  and 
of  the  surrounding  medium  is  n2  with  n2  <  ni  ,  then  the  ray  suffers 
reflections  within  the  core  of  the  guide  losing  a  fraction  of  its  power  at 
each  reflection.  An  equivalent  situation  exists  for  graded  media,  where 
a  ray  loses  power  each  time  it  touches  the  outer  caustic  or  turning  point 
along  the  ray  path  in  the  core. 

The  ray  power  transmission  coefficient  T  is  defined  as 

T  _  Power  lost  from  the  incident  ray 

Power  in  the  incident  ray  ^  J 

The  value  of  T  represents  the  fraction  of  power  lost  at  each  reflection  and 
for  an  asymmetrical  guide,  such  as  we  have  here,  it  can  vary  from  one 
reflection  to  the  next.  This  definition  is  equivalent  to  the  definition 
of  the  classical  ray  transmission  coefficients  of  Fresnel. 

The  model  for  ray  attenuation  is  formulated  as  follows.  Consider  the 
situation  where  a  ray  carrying  initial  power  PQ  at  z  =  0  undergoes  a 

fractional  loss  Ti  at  the  first  reflection,  then  after  this  reflection 
there  will  be  an  amount  of  power  Pi  left  travelling  along  the  ray,  where 

P,  =  PQ(1  -  T,)  (5) 

After  a  second  reflection  we  have  power  P2  left  in  the  ray, 

P2  =  PQ(1  -  TO  (1  -  T2)  (6) 


where  T2  denotes  the  loss  at  the  second  reflection, 
the  remaining  power  is 


After  n  reflections 


P  =  P 
n  o 


Cl  -  T.) 
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If  the  distance  to  the  first  reflection  is  z  ,  and  the  distance  between 

Pi 

the  i  -  l'th  reflection  and  the  i ' th  reflection  is  z  .  then  the  distance 
travelled  along  the  guide  by  the  ray  to  the  n'th  reflection  is 


n 


i=l 


It  is  convenient  to  express  the  power  variation  along  the  guide  in  the  form 
of 


Recasting  (9)  as 


P(z)  =  P(o)  exp  (  -7rz  ) 


-v 


=  In 


P(z) 


and  substituting  for  P(z)  and  zn  for  z  gives 


or 


'7r  )  ^i 
i=l 


In (1  -  T.) 


i=l 


zL  ln(1 " Ti) 


n 


i=l 


Only  when  T  «  1  does  (12)  reduce  to  the  simple  form 


n 


i-1 


(9) 


(10J 


(ID 


(12) 


(13) 
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The  waveguide  loss  between  the  transmitter  A  and  the  receiver  B  is 
therefore 


n 

ln(l  -  T.)J  (14) 

i»l 

if  we  chose  to  ignore  small  effects  such  as  the  Goos-Hanchen  shift(ref . 10) . 

This  then  becomes  a  powerful  method  for  solving  large  numbers  of  dielectric 
waveguide  problems. 

2.3  A  model  for  received  power 

So  far  we  have  examined  some  simple  two-dimensional  models  needed  for  a 
geometric  optics  based  solution  for  the  waveguiding  problem.  We  first  of  all 
need  to  extend  these  concepts  to  three  dimensions  before  we  can  calculate  the 
received  power  at  B. 

The  power  received  by  an  antenna  P  is  given  by 

K 


P(z)  =  P(o)exp 


I'E 


P  =  power  density  (D)  x  effective  area  (A) 
K 


where  the  aerial  capture  area  is 


A  = 


(IS) 


(16) 


Gr  is  the  gain  of  the  receiving  antenna  and  X  is  the  wavelength. 

The  power  density  within  a  flux  tube  is 

D  _  fraction  of  transmitted  power  in  flux  tube 

area  of  tube 

For  a  small  flux  tube  the  fraction  of  the  total  solid  angle  subtended  is 


(17) 


d e#> 

4ir 


and  substituting  this  in  (17)  and  using  (15)  and  (16)  we  get  the  conventional 
free  space  formula 
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Where  guidance  is  involved  then  the  cross  sectional  area  of  the  tube  is  l  Rdd> 
and 


Pq  d0/47r 


—  .  t.  R 
4JT 


where  1  is  the  width  of  the  flux  tube  and  d0  is  the  angular  width  of  the  tube 
in  the  dimension  in  which  guidance  occurs  as  illustrated  in  figure  3.  It  is 
important  to  note  that  this  implies  1  dependence  in  the  orthogonal  dimension, 

R 

ie  no  guidance  in  the  d <p  direction.  If  there  is  curvature  in  the  plane 
normal  to  the  xy  plane  defined  on  figure  10  then  additional  focussing  would 
be  achieved.  That  avenue  has  not  yet  been  investigated.  As  far  as  the  ray 
tracing  procedure  is  involved  we  see  that  it  is  necessary  only  to  calculate 
l.  If  guidance  also  occurred  in  the  orthogonal  direction  then  it  would  also 
be  necessary  to  calculate  the  depth  of  the  flux  tube  as  well. 

If  power  is  lost  from  the  flux  tube  by  either  absorption  or  tunnelling(ref . 11) 
then  the  power  received  must  be  attenuated  according  to  equation  (14)  since  it 
is  the  total  power  only  which  remains  constant  not  just  the  propagated  power. 
This  extension  of  geometric  optics  enables  us  to  take  account  of  curvature 
losses  in  the  same  way  as  Keller(ref . 9)  takes  account  of  diffraction  terms. 

We  therefore  have  a  method  of  calculating  the  received  power  within  each 
particular  flux  tube  enhanced  by  the  gains  of  the  receiver  and  transmitter 
antennae.  As  the  number  of  flux  tubes  increases  the  power  density  approaches 
its  limiting  value. 

If  the  transmitter  situated  at  A  has  the  antenna  polar  diagram  specified  in 
Appendix  IV  then  the  power  transmitted  into  the  ionosphere  can  be  described 
in  terms  of  ray  concepts,  with  the  power  associated  with  each  ray  being 
determined  from  the  antenna  polar  diagram.  In  what  follows  the  gains  of  the 
antennas  will  be  normalised  to  a  standard  dipole  (G  =G„=1.6)  but  still 

including  the  effects  of  earth  reflection  at  both  A  and  B  according  to 
Appendix  IV. 

It  is  straightforward  to  calculate  the  points  of  intersection  of  each  ray 
with  the  magnetic  field  line  and  its  locus  in  between  until  it  intersects  the 
surface  of  the  earth  in  the  opposite  hemisphere.  The  geometry  for  the 
intersection  of  a  ray  and  the  magnetic  field  line  is  written  down  in 
Appendix  V. 

At  each  point  of  intersection  the  ray  transmission  coefficient  T  is 
calculated  using  the  Generalized  Fresnel’s  Law(ref.ll)  for  reflection  from  a 
curved  dielectric  interface.  This  law  incorporates  the  geometric  optics 
extensions  needed  to  account  for  tunnelling  transmission.  Then  summing  over 
all  flux  tubes  originating  from  the  transmitter,  contours  of  the  received 
power  for  points  near  B  can  be  derived. 
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3.  RESULTS 

The  equations  described  in  the  previous  section  were  programmed  onto  a  high-speed 
digital  computer  to  calculate  the  power  received  at  B  from  a  transmitter  at  A. 

The  power  was  normalised  with  respect  to  1  kW  of  transmitter  power  and  standard 
dipoles  at  each  end  of  a  circuit  between  Darwin  and  Yamagawa  whose  coordinates 
are  listed  in  Table  1  below. 

The  powers  received  in  dB  with  respect  to  one  milliwatt  are  plotted  in  figures 
4  to  8  against  range  from  the  geomagnetic  equator  for  wavelengths  of 
X  =  3.0,  4.0,  4.5,  5.0  and  6.0  m.  The  maximum  height  of  the  magnetic  field 

line  coincident  with  the  dielectric  discontinuity  is  taken  to  be  600  km  and  the 
critical  frequency  which  is  determined  by  the  ionization  density  is  taken  to 
correspond  with  the  predictions  of  the  Ionospheric  Prediction  Service  (IPS) 

Series  U  publications  for  March  to  April  1970.  The  position  B  corresponds  to 
the  range  of  2170  km  from  the  equator. 

Experimental  results  for  this  circuit  have  been  analysed  by  Kuriki  et  al 
(refs. 12,13)  and  plotted  against  frequency.  Their  results  were  averaged  over 
two  periods  corresponding  to  the  autumnal  and  vernal  equinoxes  for  a  period  of 
several  months  during  1970.  We  estimate  the  average  critical  frequency  for 
March  to  April  as  =  10.2  Mc/s  and  for  September  to  October  as  f  =  10.8 

Mc/s.  The  results  show  an  absolute  agreement  to  within  6  dB  with  the  results 
of  our  idealised  computer  model  over  the  whole  range  of  frequencies  50  to  100 
Mc/s.  The  comparison  is  shown  plotted  in  figure  9. 

Preliminary  investigations  show  that  the  maximum  height  of  the  field  line 
required  to  give  these  results  is  about  50%  higher  than  the  height  of  the  peak 
in  the  ionization  density.  If  it  was  lower  then  the  magnitude  of  the 
dielectric  discontinuity  on  the  walls  would  be  insufficient  to  maintain  guidance 
while  if  it  was  higher  then  the  ionization  at  the  top  of  the  field  line  would  be 
insufficient. 

As  yet  the  reason  for  the  bulge  in  the  experimentally  observed  results  near 
88  MHz  has  remained  obscure  but  this  bulge  was  not  observed  in  the  results  of 
Bowen  et  al(ref. 14)  or  Nielson(ref . 1)  in  their  results.  However,  the  variation 
of  received  power  with  frequency  is  not  linear  as  proposed  by  Kuriki  et  al(ref. 13) 
even  though  over  a  limited  range  of  frequencies  it  may  appear  to  be.  He  used  a 
theory  based  on  scattering  from  ionospheric  irregularities  but  needed  to  make 
questionable  assumptions  regarding  the  nature  of  scattering  from  these 
irregularities  and  their  frequency  of  occurrence.  His  scattering  theory 
produced  the  straight  line  with  f^10  dependence  superimposed  on  figure  9. 

In  contrast  the  calculations  used  in  this  model  attribute  the  power  to  guidance 
rather  than  scattering  and  rely  merely  on  the  existence  of  field  aligned 
irregularities  that  are  likely  to  be  associated  with  rising  plasma  bubbles. 
Discrepancies  between  observed  and  calculated  results  could  easily  be  due  to 
imperfections  in  the  shape  or  continuity  of  the  elliptical  surface  assumed  in 
this  model.  Alternatively,  if  the  ionization  within  the  bubble  was  not  totally 
depleted  then  a  slightly  smaller  critical  angle  would  be  associated  with  points 
on  the  ellipse  and  this  would  reduce  the  received  power. 

Changes  in  the  critical  frequency  of  6%  observed  between  the  autumnal  and  vernal 
equinoxes  produces  only  a  ^  dB  change  in  the  received  power  for  the  Darwin/ 
Yamagawa  circuit  at  X  =  6m  using  this  geometry  rather  than  the  3  dB  observed 
experimentally. 


I 
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Since  the  results  of  this  initial  investigation  are  so  close  to  the  experimental 
results  further  analysis  appears  to  be  warranted.  In  a  subsequent  study  we  will 
perturb  all  the  parameters  associated  with  this  model  to  determine  their  allowable 
range  and  to  seek  an  experimental  method  of  confirming  their  behaviour. 


4.  CONCLUSIONS 

A  simple  model  for  the  received  power  on  a  transequatorial  circuit  described  by 
figure  1  has  been  proposed  which  may  explain  night-time  radio  wave  propagation  in 
the  ionosphere.  The  mechanism  has  been  ascribed  to  guidance  along  field  aligned 
irregularities  associated  with  rising  plasma  bubbles  in  the  ionosphere.  Good 
agreement  has  been  demonstrated  between  experimental  results  obtained  on  a  circuit 
between  Darwin  and  Yamagawa  and  theoretical  results  for  an  ideally  formed  guiding 
irregularity. 

Further  investigation  is  required  to  optimise  the  model,  to  determine  an 
experimental  technique  of  validating  it  and  to  evaluate  the  effect  of  imperfec¬ 
tions  in  the  guiding  surface. 
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NOTATION 


intercept  of  ray  with  y  axis 
variable  defined  in  Appendix  V 
critical  frequency 
plasma  frequency 

maximum  height  of  field  line  chosen 


wavenumber 


component  of  wave  vector  I  kl 
ray  slope 

incident  ray  slope 

slope  of  ellipse  at  point  of  reflection 
reflected  ray  slope 
refractive  indices 

polar  coordinate  for  point  on  field  line 

minor  semi  axis  of  ellipse 

major  serai  axis  of  ellipse 

ionospheric  scale  height 

antenna  height 

electron  density 


radius  of  earth 


ray  transmission  coefficient 


Fresnel  coefficient 


coefficient  of  «  Chapman  layer 
dielectric  permittivity 
wavelength  of  source 
critical  angle 
angle  of  incidence 

radius  of  curvature  of  ellipse  at  point  of 
reflection 
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phase  path  difference 

ray  direction  from  transmitter 

variable  defined  in  Appendix  IV 
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APPENDIX  I 

STRUCTURE  OF  THE  EQUATORIAL  FIELD 

The  equation  of  a  metnetic  field  line  from  a  point  magnetic  dipole  located  at  the 
centre  of  the  earth  is  given  by  the  equation (ref . 15) 

T 

2a  =  constant  (1.1) 

COS  V 

where  r,  0  are  the  polar  coordinates  of  a  point  on  the  field  line  as  shown  in 
figure  10.  The  intersection  of  a  field  line,  which  passes  through  a  maximum 
height  h  at  the  geomagnetic  equator  (G) ,  with  the  surface  of  the  earth  (radius 
R)  occurs  where  latitude  0  is  given  by 


6  =  ±arccos  (1.2) 

In  the  local  plane  cartesian  coordinates  Oxy  position  of  the  intersection  is  at 
(Rsin0 ,  Rcos0  -  R/2). 

If  the  discontinuity  is  assumed  to  begin  at  conjugate  points  on  opposite  sides 
of  the  equator  S,  S'  at  a  height  of  100  km  then  in  local  cartesians  Oxy  the 
points  S,  S'  using  equation  (1.2)  are  at  (R  +  100)  sin0 ' ,  (R  +  100)  cos0 ' ,  - 
R/2,  where  0'  is  the  latitude  of  intersection  of  the  field  line  contour  and  the 
contour  100  km  above  the  earth's  surface.  Fitting  an  ellipse  to  the  points  G, 
S,  S'  gives  a  best  fit  curve  having  the  equation 


_ X _  ") 

(R  +  100)  sin0^. 

^  j 

J  i  .  /( R  +  1OO)cos0  '  -  R/2  \  2'J 

V  \ - (TTh)72 - / 

♦  [  _ L _ 

V  (R  +  h)/2 


1 

(1.3) 


In  the  region  of  the  equator  the  error  involved  in  assuming  the  equation  of  the 
field  line  to  be  equation  (1.3)  for  h  =  500  km  is  less  than  1%  compared  to 

equation  (I . 1)  over  the  whole  of  the  curve  above  the  ground. 
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APPENDIX  II 

MODEL  FOR  ELECTRON  DENSITY 

An  appropriate  electron  density  model  is  the  Chapman  model  defined  by  = 
Reference  16.  The  frequency  f  is  defined  in  terms  of  the  critical  frequency 

fc  by 

fN  =  fc  6XP  (<X(1  "  Z  '  6"Z)) 


where 


z 


h  -  h  max 
H 


a  =  0.5  for  an  a  Chapman  layer 


h  max  is  the  height  of  the  electron  density  peak. 

A  plot  of  height  against  ionization  density  for  varying  values  of  H  with 
f  =  12  MHz,  h  max  =  400  km  appears  as  Figure  11.  We  choose  H  =  60  km  to 

give  the  distribution  of  Booker (ref . 2) . 

The  electron  number  density  N  is  then  given  by  the  formula 

N  =  1.24  x  10"*  f?t 

N 


J 
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APPENDIX  III 

CHANGE  IN  REFRACTIVE  INDEX  AT  AN  IRREGULARITY 

The  fractional  change  Ac  in  the  dielectric  permittivity  e  in  an  ionized  medium  is 
related  to  the  change  in  ionization  density  An  by  the  relation 


_1_  2 

w2  * 
it  e 


where  r  is  the  classical  radius  of  the  electron 
e 


(2.82  x  10_ls  m) 

Strictly,  X  should  be  the  wavelength  taking  into  account  the  ionization  density 
N,  but  for  VHF  propagation  in  this  model  we  will  take  X  to  be  the  free  space 
wavelength.  At  the  irregularity  therefore  we  will  see  a  change  in  the  magnitude 
of  the  refractive  index  discontinuity  with  height  and  a  corresponding  change  in 
the  critical  angle  0  where 


6  =  arcos  — - 

c  n2 


where  nj ,  m  are  the  refractive  indices  on  either  side  of  the  discontinuity. 

In  the  case  of  a  totally  depleted  plasma  bubble  rising  in  the  ionosphere  we  put 
ni  =1. 

If  we  arrange  =  1  then  n2  is  determined  and  hence  0  . 

N  C 
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APPENDIX  IV 

ANTENNA  POLAR  DIAGRAMS 

The  polar  patterns  characterising  either  a  tiansmitting  or  receiving  antenna  at 
height  Hi  above  the  ground  are  assumed  to  be  identifical .  The  received  field 
is  given  by  the  sum  of  the  direct  and  reflected  wave  amplitudes  as  depicted  in 
Figure  12.  At  low  angles  we  can  assume  that  the  reflection  coefficient 
of  the  surrounding  earth  is  -1  and  if  the  wave  amplitude  is  A  then  the  total 
amplitude  received  assuming  unity  antenna  gain  is 


=  Asinwt  -  A  sin(wt  +  <t> ) 


(IV. 1) 


where 


0 


(IV. 2) 


and 


X  =  operating  wavelength 


This  reduces  to 


A  =  angle  of  propagation. 


Ajsinwt(l  -  cos0)  -  coswt  sin0j 


=  A'  sin(wt  -  4>) 


where 


<t>  = 


(IV. 3) 


The  magnitude  of  the  power  from  a  transmitting  antenna  in  direction  A  is  then 


A'2  = 


A2  (2  -  2cos0) 


(IV. 4) 
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APPENDIX  V 


THE  GEOMETRY  OF  REFLECTION  FROM  AN  ELLIPTICAL  SURFACE 

Simple  geometry  as  shown  in  figure  13  is  used  to  determine  the  angles  of 
incidence  and  reflection  at  a  dielectric  discontinuity.  If  the  equation  of  the 
ellipse  is 


A2 


O' .  -  YI2 

bt 


i 


and  of  the  incident  ray  is  y  =  mx  +  b  where  b  is  determined  by  the  initial 
conditions  then  the  point  of  intersection  is  given  by  the  solution 

-A2  me  +  AB'/B2  -  c2  +  A2m2 

x  =  - (FT??) 


where 


c  =  b  -  Y 


The  slope  of  the  ellipse  at  the  point  of  intersection  is  m2  given  by 


-xB 

1,12  =  AV  -  x2 

If  the  angle  of  incidence  measured  with  respect  to  the  reflecting  surface  is 
0 i  *  arctan  m2  -  arctan  mi  and  the  slope  of  the  reflected  ray  is  m3  then 

m3  *  tan  0 

r 

and 


m3  -  2v\2  ~  nil 


* 

k 
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TABLE  1.  COORDINATES  OF  DARWIN  AND  YAMAGAWA 


GEOGRAPHIC  COORDINATES 

LATITUDE  RELATIVE  TO 

GEOMAGNETIC  EQUATOR 

DARWIN 

130°  50' 

E,  12°  20'  S 

19°  32'  S 

YAMAGAWA 

130°  37' 

06"  E,  30°  12'  05"  N 

23°  O'  N 

UOO  1600  1800  2000  2200  2400  2600  2800 

RANGE  (  km) 

Figure  6.  Plot  of  normalised  received  signal  strength  against  distance 
from  the  equator  for  X  =  4.5  metres  and  h  =  600  kilometres 
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Figure  7.  Plot  of  normalised  received  signal  strength  against  distance 

from  the  equator  for  X  =  5.0  metres  and  h  =  600  kilometres 


km) 

Figure  8.  Plot  of  normalised  received  signal  strength  against  distance 
from  the  equator  for  X  =  6.0  metres  and  h  =  600  kilometres 
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Figure  9 


Figure  9 


Plot  of  normalised  received  signal  strength  against  frequency 
for  both  theoretical  and  experimental  results  as  discussed  in 
the  text 
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Figure  11.  Plot  showing  relat: 
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